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(57) Abstract 

In order to optimize the determination of radio-frequency energization in an NMR experiment, the Schrodinger 
equation is written relative to the amplitude of probability of the wave function of spins resonating at a given resonance 
frequency and of which the Hamiltonian depends on the frequency and on the radio-frequency energization to which they 
are subjected. Analytical polynomial, expressions are deduced therefrom for the probability amplitudes at the end of the 
energization as a function of the resonance frequency. The coefficients are optimized by means of calculation algorithms, 
from which the characteristics of the energization are deduced. 

(57) Abrege 

Pour permettre une optimisation de la determination d'une excitation radiofrequence dans une experimentation de 
RMN, on ecrit l'equation de Schrddinger relative a l'amplitude de probabilite de la fonction d'onde de spins resonnant a 
une frequence de resonance donnee et dont l'hamiltonien depend de la frequence et de 1' excitation radiofrequence qu'ils 
subissent. On en deduit des expressions analytiques polynomials des amplitudes de probabilite a la fin de 1'excitation en 
fonction de la frequence de resonance. On optimise les coefficients a I'aide d'algorithmes de calcul et on en deduit les ca- 
racteristiques de 1'excitation. 
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PROCEDE D'EXCITATION RADIOFREQUENCE 
DANS UNE EXPERIMENTATION DE RMN. 



La presente invention a pour objet un procede d'excitation 
radiofrequence dans une experimentation de resonance magnetique 
nucieaire (RMN). Elle trouve plus particuiierement son application 
dans ie domaine medical ou des machines de RMN sont utiiisees a 
des fins de diagnostic chez l'homme. Elle peut neanmoins trouver 
son application dans les domaines du laboratoire ou de 1'industrie.des 
que des phenomenes de resonance magnetique sont mis en jeu. 

Le phenomene de la resonance magnetique resulte de l'orien- 
tation prise par les moments magnetiques de particules cfun corps 
quand ce corps est soumis a un champ magnetique orientateur 
constant. Pour faire apparaitre ce phenomene on excite le corps par 
une excitation radiofrequence. On peut mesurer le signal de des- 
excitation des particules quand elles retournent a leur equilibre, 
quand leur moment magnetique se reoriente avec le champ orien- 
tateur, a la fin du basculement de cette orientation dependant de 
l'excitation. L'excitation est appiiquee a la totalite du corps soumis 
a examen, et toutes les particules du corps emettent un signal de 
desexcitation a son issue. Dans le domaine medical ou on cherche a 
representer des images de coupes du corps examine (d'un patient), on 
a coutume de cantonner la naissance du phenomene de resonance a 
des tranches selectionnees dont on veut donner l'image. Dans ce but 
on modifie localement les conditions de resonance de maniere a ce 
que seule la tranche selectionnee se trouve dans des conditions de 
resonance correspondant a la frequence de l'excitation radio- 
frequence. En dehors de la tranche, les moments magnetiques ne 
sont pas bascules, les particules exterieures a la tranche ne res- 
tituent aucun signal de desexcitation. 

La modification des conditions locales de resonance s f obtient 
en augmentant, ou en diminuant, Tintensite du champ orientateur 
dans la tranche selectionnee. Pour des raisons scientifiques et 
technologiques il n ! est pas possible de modifier brutalement 
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l'intensite cfun champ magnetique continu de part et tfautre rfune 
tranche de l'espace. La modification du champ prend alors la forme 
d'une variation reguliere en fonction de 1'abscisse des iieux de 
l'espace sur un axe perpendiculaire a la tranche a selectionner. 

5 L'intensite du champ orientateur subit done selon cet axe un 
gradient, et pratiquement les champs supplementaires, additionnes 
au champ constant et homogene dans tout l'espace pour provoquer 
cette variation, sont commune ment appelles ies gradients de champ. 
La frequence de resonance du phenomene de resonance est propor- 

lO tionnelle a Tintensite du champ orientateur. Pour selectionner une 
tranche il convient done, dans 1'excitation radiofrequence, que le 
spectre de 1'excitation ne contienne des composantes que comprises 
dans une bande limitee par des frequences de coupures corres- 
pondant aux frequences de resonance des particules des bords de la 
tranche. 

L'obtention cfun signal d'excitation a spectre etroit entraine 
que ce signal ne peut etre instantane. II dure un certain temps. En 
pratique des puissances de quelques kilowatts, appliquees pendant 
des durees de quelques dizaines de millisecondes, sont courantes. Le 

20 fait que Timpulsion cfexcitation n ! est pas instantanee et qu ! au 
contraire elle dure un certain temps provoque alors des pertur- 
bations dans le signal de desexcitation. En eff et, du fait du gradient 
de champ, les particules des bords de la tranche ne resonnent pas a 
des memes frequences. Par exemple si le milieu de la tranche 

25 resonne a une frequence f Q , les bords amont et aval de la tranche, 
par rapport a Torientation du champ, peuvent resonner a des 
frequences f q - A f . ou f q + A f respectivement et il en resulte que 
i'historique de Texcitation dans les bords de la tranche rfest pas le 
meme qu'au centre bien que la duree de fexcitation ait ete la meme 

30 pour toutes les particules de la tranche. Si le basculement geo- 
metrique de tous les moments magnetiques est correctement obtenu 
pour une energie cfexcitation donnee (par exemple lis se retrouvent 
tous dans un plan perpendiculaire au champ orientateur pour une 
excitation a 90°), leur orientation rfest alors pas unique. Elle peut 
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au contraire etre distribute dans toutes les directions contenues 
dans un plan. L'orientation de chaque moment y depend de Tabcisse 
des points consideres dans la tranche. En effet dans leur mouvement 
de precession forcee pour rejoindre le basculement geometrique 
5 recherche, les moments magnetiques ont emprunte des trajectoires 
differentes. Ces trajectoires sont differentes car les vitesses de 
precession sont liees aux frequences de resonance. Comme i'orien- 
tation initiale des moments magnetiques est la meme pour tous, on 
comprend qu'il n'y a aucune raison a priori pour que cette orien- 

10 tation soit unique a la fin du basculement. Si on n'y prend garde les 
moments magnetiques des divers points de la tranche peuvent se 
trouver a ce point desorientes, dephases, les uns par rapport aux 
autres que le signal electromagnetique de restitution peut avoir une 
moyenne nuile. II est alors impossible d'en tirer une information 

15 physique. 

En outre le phenomene de resonance magnetique ne se situe 
pas a une frequence unique, mais plutot a des frequences voisines. 
L'ecart entre ces frequences et une frequence de reference est 
appele le decalage chimique (en literature anglo-saxonne : chemical 

20 shift). Cet ecart est en general inferieur a 2 A f. II est recherche 
par les praticiens une excitation qui serait capable, pour une 
experimentation donnee, de differencier la reponse correspondant a 
chaque valeur de se decalage chimique. De cette maniere il serait 
possible, en particulier dans les images medicales, de faire une 

23 separation eau-graisse dans les tissus. En effet pour les differents 
tissus, si le decalage de I'eau par rapport a la graisse a une meme 
valeur en frequence, les amplitudes des composantes n ! ont pas des 
memes valeurs d*un tissus a Tautre. Par ailieurs, dans une autre type 
^application et en particulier dans le domaine de Timagerie, il est 

30 envisage de construire des excitations dites multicoupes ou plutot 
multi tranches. Le but serait de profiter d'un gradient qui s'etend 
dans i'espace sur une longueur bien superieure a Tepaisseur de la 
tranche et de disposer d'une excitation susceptible rfexciter, en 
meme temps, plusieurs tranches, non contigues de preference, de 
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maniere a acceierer les procedes de reconstruction cFimage. Or le 
travail entrepris pour la definition des impulsions d'excitation 
classique (a 90° et a 180°) a paru a ce point complexe que i'homme 
du metier n'envisage pas rfen augmenter encore la complexite pour 
produire maintenant de telles excitations. 

Dans Tetat de la technique le probleme de la determination du 
signal d'excitation a ete resolu par des approximations successives 
et par des tatonnements qui peu a peu ont permis de determiner 
Thistorique de son amplitude pendant la duree de Timpulsion de 
maniere a ce que des resultats desires soient atteints. Cependant, 
malgre les mois de travail necessaires pour arriver a ces resultats, 
Toptimisation des impulsions d'excitation n'est pas garantie, 
notamment en ce qui concerne la puissance totale mise en oeuvre 
ainsi que la puissance instantanee. Le but de la presente invention 
est de proposer un procede de determination de Texcitation qui 
permette en plus de tenir compte d'un critere cPoptimisation, meme 
eventuellement different de ceux-ci. 

La technique utilisee dans Petat de la technique pour 
determiner les excitations consiste essentiellement a diviser une 
2Q duree predeterminee (Jexcitation en un nombre suffisant de 
periodes, par exempie ^0, et a calculer pour chacune des abcisses de 
la tranche la rotation, pendant chaque periode, du moment magne- 
tique des protons qui s'y trouvent. Le calcul est entrepris en 
fonction cfun parametre representatif de Tamplitude de i'excitation 
25 pendant la periode correspondante. On fait evoluer dans le temps ce 
parametre pendant la duree de Texcitation. Le mouvement du 
moment magnetique est un mouvement de precession du a la 
frequence de resonance, done dependant de l'abscisse dans la 
tranche, combine a un mouvement de nutation impose par l'exci- 
2Q tat ion radiofrequence, identique pour toute la tranche, Aussi il met 
en jeu pour determiner la composante selon une direction de 
Tespace du moment magnetique cfune particuie a la fin cfune 
periode, les composantes selon les trois directions de l'espace du 
moment magnetique de cette particuie avant le debut de la periode, 



ces trois composantes dependant, du fait du gradient, de i'abcisse du 
point dans Ja tranche. En consequence le calcuJ doit etre entrepris 
abcisse par abcisse pour toutes les abcisses dans i'epaisseur ie long 
de la tranche. 

Dans la pratique ces calculs ont ete entrepris pour une 
centaine d'abcisses reparties regulierement au trayers, et un peu de 
part et cPautre, de la tranche. Le travail de determination a 
consiste a faire la totaiite de ce calcul, et a le reprendre un certain 
nombre de fois en modifiant, periode par periode, Thistorique de 
revolution du parametre de maniere a obtenir en fin compulsion une 
definition recherchee en amplitude et en phase du basculement des 
moments magnetiques des particules soumises a ^experimentation. 
Ces calculs concernent toutes les particules dans la tranche, et 
voisine de la tranche, c'est a dire un ensemble de frequences de 
resonance. Dans une autre approche, justifiee uniquement pour de 
faibles angles de nutation, par exemple inferieurs a 30°, on utilise 
le fait que l'expression de l'aimantatipn transverse en fonction de la 
frequence de resonance est directement la transf ormee de Fourier 
du signal temporel d'excitation. Cette hypothese devient peu valide 

■ 

quand l f angle de nutation recherche vaut 90°. Elle est 
completement fausse quand il vaut 180°. 

La presente invention permet de remedier aux inconvenients 
cites en proposant une methode de calcul de 1'impulsion d'excitation 
completement differente. Dans cette methode on ne calcule pas en 
chaque abcisse directement V angle de nutation de chaque moment 
magnetique. On exprime seulement, sous une forme analytique poly- 
nomiale en fonction de la frequence de resonance, la probability que 
ce moment ait une orientation desiree a la fin de 1'excitation. 
L'expression sous forme poiynomiale permet de s ! aff ranchir, dans les 
calculs, de Tinformation relative a la localisation des protons 
concernes (ou de leur frequence de resonance ce qui revient au 
meme) de maniere a pouvoir proposer une solution presentant pour 
toutes les frequences de resonance impliquees un meilleur com- 
promis vis a vis d'un critere ^optimisation. En ramenant cette 
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optimisation a des syntheses de filtres par transformee de Fourier, 
Tinvention permet alors d'aboutir rapidement a une solution re- 
cherchee (en quelques minutes). 

La presente invention a pour objet un procede cfexcitation 

5 radiofrequence dans une experimentation de RMN, excitation en 
reponse a laquelle les bascuiements de moments magnetiques de 
particules soumises a l'exper indentation sont definis en amplitude et 
en phase pour un ensemble de frequences de resonance donne 
caracterise en ce que : 

10 - on ecrit l'equation de Schrodinger relative a 1'amplitude de 

probalite de la f onction rfonde de spins resonnant a une frequence de 
resonance donnee et dont Thamiltonien depend de la frequence et de 
Pexcitation radiofrequence qu f ils subissent, 

- on integre cette equation par rapport au temps, 

13 - on en deduit les expressions analytiques polynomiaies des 

amplitudes de probabilite a ia fin de i'excitation en fonction de la 
frequence de resonance, ces expressions polynomiaies etant de degre 
egal a un nombre de periodes elementaires dont la somme des durees 
est egaie a la duree de Texcitationi, et ces expressions etant - 

20 parametrees par des coefficients representatifs de Thistorique de 
^excitation pendant ces periodes, 

- on optimise ces coefficients a Taide d'aigorithmes de calcul 
de fa^on a ce que les polynomes ^approchent au mieux d ! une forme 
desiree 

25 - on deduit les caracterisitiques, pendant chaque periode, de 

l'excitation qui conduit a cette reponse definie. 

- et on excite les particules avec une telle excitation. 

La presente invention sera mieux comprise a la lecture de la 
description qui suit et a Texamen des figures qui i'accompagnent. 
30 Cette description n'est pas limitative de Tinvention. Notamment, 
^application cfimagerie evoquee n'est qu'un cas particulier de de- 
termination d f une excitation radiofrequence. 

Les figures montrent : 

- figure 1 : une representation schematique tfune machine de 
RMN ; 
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- figures 2 et 3 : des allures cfenveioppes cfampiitude 
d'excitation determinees selon Tinvention. 

La figure 1 represente schematiquernent une machine 
d 1 imager ie par resonance magnetique nucleaire. Cette machine 

5 comporte essentiellement des moyens 1 pour produire un champ 
magnetique orientateur, constant et homogene, B Q dans un lieu ou 
est place un corps 2 a examiner. Le corps 2 est soumis a cet endroit 
a une excitation radiofrequence produite par generateur 3 en 
relation avec une antenne 4. L'antenne sert egalement a prelever le 

10 signal de desexcitation et a Tacheminer, au moyen d'un duplexeur 5, 
sur des moyens de reception et de traitement 6. A i'issue du 
traitement, l'image d f une coupe 7 du corps 2 peut etre visualisee sur 
un dispositif 8 de visualisation. Pour soumettre les moments magne- 
tiques des protons du corps 2 situes dans la tranche 7 au phenomene 

15 de resonance, et eux seuls, on a mis en oeuvre des bobines de 
gradients ? alimentees par un generateur 10 d'impulsions de 
gradient. Tous ces moyens fonctionnent sous ie pilotage cfun 
sequenceur ii. En supposant que !e champ magnetique suple- 
mentaire applique par les bobines 9 varie line aire ment selon un axe 

20 2 paralleie au champ Bq, on peut dire que les bords de la tranche 7, 
bordsi2 et 13 respectivement en amont et en aval, rencontrent des 
conditions de resonance differentes et que les frequences de 
resonance des protons qui s ! y trouvent sont differentes. 

Dans 1'invention, pour determiner Texcitation radiofrequence, 

25 on ecrit ^equation de Schrodinger relative a Tamplitude de proba- 
biiite de la fonction cfonde de spins resonnant a une frequence de 
resonance donnee. Cette expression est la suivante : 



30 



Dans cette formule H est Phamiltonien de I'equation, il 

a 

depend de la frequence et de Texcitation radiofrequence que 
subissent les spins re sonants. Le terme | ^> represente la fonction 
d'onde. Cet hamiltonien peut s'ecrire de la maniere suivante : 
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-w -w + iw 

x y 

- w - iw w 
x ' y 

2 

Dans ces equations j est le nombre complexe tei que j = -1. Les 
pulsations w, et w et w sont, respectivement, fecart de pulsation 

* y 

de resonance dans la tranche 7 en fonction de Tabcisse selon z des 
protons concernes, et les representations des composantes radio- 
frequence du champ d 1 excitation emis par Tantenne *. Ces 
composantes sont evaluees dans un repere tournant a une pulsation 
Wq, correspondant a la porteuse de I'excitation radiof requence et a 
la frequence de resonance des protons situes au centre de la tranche 
7. De ^equation de Schrodinger on tire Texpression suivante : 



a~ 2 



15 



20 



25 



30 



da 




jw 


dt 




~2 


db 






dtl 

1 1 


i ' 


2 

— 



w 1 e 



-j<t> 



-JW 



I 



Ou a et b sont les amplitudes de probability pour que les spins des 
protons consideres soient dans un etat respectivement positif ou 
negatif apres transformation par i'operateur moment cinetique selon 
z. Dans l'expression I on a pose : 



j (w + jw ) = w. e ' ^ 
J • x ; y 1 



ou w^ est I'amplitude de I'excitation radiof requence et ou <|) est sa 
phase cfe mission par rapport aux axes x et Cette phase est 
conventionellement nulle quand remission est selon l'axe -y. 

Pour integrer l'equation I par rapport au temps on procede 
d ! une maniere preferee a une approximation, de type connu, dite par 
"impulsions dures" ou encore "hard pulse". Selon cette approximation 
on considere que les protons sont soumis, pendant I'excitation, a 
deux actions differentes. Une premiere action permanente est 
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15 



25 



relative a la presence du gradient de champ orientateur dans la 
tranche 7. Elle a pour effet le decalage de pulsation de resonance w. 
La deuxieme action correspond a Texcitation radiofrequence 
proprement dite. Dans la pratique evoiue dans le temps et c'est 
son evolution que Ton recherche. En echantillonnant on peut la 
representer sous la forme rfune suite de valeurs bloquees en paliers 
pendant des periodes elementaires et evoluant brutalement, d'un 
palier a un autre, de la fin June periode au debut d'une autre, 
L' approximation evoquee consiste a remplacer revolution par palier 
de la deuxieme action par une succession compulsions de type Dirac, 
de duree i tres courte et dont Tamplitude w. est egale a (t), x/6 
ou (t) est la valeur des paliers de l'excitation radiofrequence et 
ou t est la duree cfune periode elementaire. L' effet de cette 
approximation est de simplifier notablement Pequation L En effet on 
peut dire que pendant les impulsions dures, pendant les diracs, elle 
s'ecrit : 





da 




• 


dt 




20 


db 






dt 





w. n 
i e 



w. 
2 



-it, 









b 



parce que w. est tres superieure a l'effet du gradient de champ 
orientateur. 

Hors des periodes d'application des excitations en Dirac, cette 
expression peut s'ecrire simplement: 



30 



da 
dt 



2 



db 
dt 



2 



puisque dans ce cas I'impulsion de Dirac n'est pas presente. On peut 
integrer chacune des deux expressions precedentes par rapport au 
temps. Le calcui est classique, et on peut ecrire respectivement : 
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Dans ces expressions ies amplitudes de probability sont affectees du 
signe - ou du signe + selon qu'elles correspondent a ces probabilites 
1^ avant ou apres 1'appli cation du dirac _i correspqndant. EUes sont 
affectees d'indices i ou i + 1 selon ie rang ± de i'impulsion de Dirac 
qui les precede. Dans ces expressions 9j = w..£ = (t).x. Si par 
ailieurs on pose pour simpiifier les expressions : 



20 



cos O.J2 



c. 
i 



et 



25 



sin 6./2 



s. 
i 



on peut ecrire qu'apres une n ieme impulsions de Dirac, en omettant 

■ 

les signes + : 



30 



n 



n 



c - s 
n n 



s c 
n n 



z 112 0 



0 z 



-1/2 



l n-l 



n-1 



De cette ecriture on va deduire des expressions de a^ et b^ en 

fonction de a^. - et b^ _ ^. Et dans ces dernieres expressions, en 

remplacant a et b par : 

n n r 
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a = 



a / z 
n 



n/2 



et 



b /z 
n 



n/2 



on peut ecrire : 



a = c a - s z B , 

n n n - 1 n / n - 1 



- 1 * 

ft = s a ,+cz B . 
rn n n - 1 n fn - 1 



10 



2 



Si on suppose connue fomentation du moment magnetique avant une 
premiere impulsion de Dirac, en particulier s'il est oriente paral- 
lelement a B Q , on peut dire que : 
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Dans la derniere expression a Q et p Q dependent de c Q et s Q qui sont 
representatifs de l'amplitude du premier dirac numerote 0. Cette 
formulation peut encore s'ecrire : 



et 



u 0 -0 
°0 = °0 = h 0 2 



0-0 

h = s o = g o 2 



oil z° vaut conventionnellement 1. Apres la deuxieme impulsion 
(numerotee n = 1) on peut ecrire : 
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a i = c i c o - ~ s i s o z ' 1 =h o z "° + h i 2-1 
h = s i c o " c i s o 2-1 = So 2 ~° + «} 2-1 

De meme apres une troisieme impulsion (numerotee n = 2) on peut 
5 ecrire : 

a 2 = c 2 ( Cl c Q - Sj s Q 2* *) -l 2 z- 1 (s 1 c Q - c t s Q z" l ) 



10 



15 



20 



soit a 2 = z" 0 + h* z 1 + h 2 z " 2 

et Ton pourrait ecrire de meme pour jS 2 : 

a « 2 ,"° 2 -1 2 -2 
£ 2 = g Q z + %x z + g 2 z 

Ainsi de suite et plus generalement apres la n plus unieme impulsion, 

numerotee n, a et B s'ecrivent i 

~ n * n 



„ _ h n -0 . n - 1 n -n 

a = n z + h. z + ... + h z 
n 0 1 n 



a n -0 n -1 n -n V 

fn = h 2 + «1 2 ■ + - + §n 2 



Avec cette expression on peut faire deux remarques 
importantes. Premierement aet^ representent, a un changement de 
variable pres, les amplitudes de probability de la fonction cfonde 
transformee par i'operateur moment cinetique selon l'axe z. On peut 

25 cfailleurs montrer, qu'au meme titre que le carre du module de a 
additionne au carre du module de b vaut 1, que le carre du module 
de additionne au carre du module de ji n vaut 1 egalement. En 
effet ia probability pour que le moment magnetique existe vaut 
toujours 1. Deuxiemement on remarque que Texpression de ct^ et p 

30 est une expression polynomiale analytique en z. Or z correspond a la 
pulsation de resonance en fonction du lieu (a cause du gradient) des 
protons concernes dans Tespace excite. Done z est dans l'expression 
du moment magnetique, le terme representatif de la frequence de 
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resonance. Au meme titre que c Q et s Q etaient determines par 
1'amplitude du premier dirac ^excitation on peut dire que ies 
coefficients h" et g. sont determines par ia succession des differents 
diracs, jusqu'au dirac numero n . Ces coefficients sont 
5 cependant exprimes, eux, independemment du lieu des protons 
concernes c f est a dire independemment de la frequence. Autrement 
dit ia connaissance de tous les h 1 ? et g" per met de deduire failure du 
moment magnetique a la fin de V excitation quelle que soit la 
frequence de resonance, c'est a dire quelle que soit la place des 

1 0 protons dans la tranche. 

Dans ^invention on se preoccupe aiors de rechercher de teis 
coefficients rfj et g^ en imposant des contraintes lies aux valeurs 
de w. En queique sorte on va imposer que les expressions de et 
p^soient les equations d'un filtre ayant une valeur donnee dans sa 
15 bande passante et ayant une autre valeur donnee a l'exterieur. 

On peut ecrire en reportant la valeur _ . et de £ R _ ^ en 
fonction des coefficients " 1 et g 1 ? " 1 de i'expression V dans 
Texpression IV de et p^, eux-memes ecrits en fonction des 
coefficients et g? : 

20 

n ,n - 1 - n - 1 

h i = V h i - S n Si . i 

et g n = s .h^ 1 + c , g n ~; 

5 i n i n °i - 1 

25 L'expression IV peut en outre etre inversee et devenir : 



30 



v 



a , = c a + s fi 
n-i n n nrn 

fi t = z(-s a + c 8 ) 
"n-1 n n n/n 



V 



A ce stade du calcul mathematique on peut faire les consi- 
derations suivantes. Si on s'est donne a n et p n , representatifs <f un 
moment magnetique a la fin de Texcitation, on peut connaissant c n 
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et representatifs du dernier dirac de l'excitation, deduire les 

valeurs du moment magnetique, a ^ et P n _^ juste avant i'appli- 

cation de ce dernier dirac. De proche en proche, en determinant a 

chaque fois le dirac precedent, on peut remonter revolution du 

moment et determiner, a l'envers, Thistorique de l f excitation depuis 

son origine jusqu'a sa fin. La question est aiors de savoir si on peut 

choisir c ou s comme on le yeut ou au contraire si ces valeurs sont 
n n 

imposees. Que ces valeurs soient imposees revient a dire que le 
dirac numero n est connu en amplitude. En examinant VIII et en se 
rappeliant que et sont des polynomes de degre n en z * il vient 
im mediate ment : 



15 



20 



et 



0 = c h n + s g n 
n n n °n 

,n n 
0 = ~ s n V c n g 0 

- 1 



En eff et le polynome de degre n - 1 en z" representatif de c* n ^ ne 

peut contenir des termes de degre n en z . De meme le polynome 

-1 

de degre n-1 en z representatif de ^ ne peut comporter des 
terr: es en z*. Ce couple de conditions revient a poser les equations 



suivantes : 



25 



30 



i e^ ^n. sin 9 L 

n n u 

r cos & 12 

n n 



n 



2 



g 



n 
n 

n 
n 



g 



n 
0 

n 



x 



Ceci revient a imposer que G et 6 soient les deux inconnues d*un 

n n 

syteme de * equations a deux inconnues. En principe on devrait se 
heurter igi a une impossibility de resolution. 

Mais le fait de dire que la probability que le moment 
magnetique existe vaut 1 revient a ecrire que : 





2 




2 _ 


a 


+ 


h 




n 
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ou encore : 



a a 
n n 



fnPn = 



1 



2 



Cette expression VII peut etre developpee en remplagant <* n et £ n 
par leur valeur donnee dans l'expression V« En effectuant ce dernier 
calcul et en classant les coefficients obtenus selon les puissances de 
z auxquelles Us s'appliquent on peut ecrire : 



u n un n n 

h 0 h n + §0 g n 

10 7n n 7n n n n n n 

h 0 h n-l + h l h n + g 0 §n-l + §1 §n 



= 0 termes en z 



-n 



= 0 termes en z 



= 0 



-(n- 1) 



15 



h Q h Q + hj hj + 



,+ h n + 
n n 



25 



30 



n n n n 
W U W U -i -1 



20 _ — 



"0 n 



n 



n n 



« > » 



• ••• 



= 1 termes en z 



= 0 



-0 



- 0 termes en z 



n 



Ce systeme cFequations comporte n+1 equations, et non pas 2n 
equations, a 2n + 2 inconnues. En effet les equations dont les 
termes sont en z l ou en z 1 sont les memes. Par contre il y a bien 2n 
+ 2 termes en h 1 ? et g. qui sont inconnus. On remarque que la 
premiere et la derniere equations reduisent Timpossibilite de deter- 
mination de 9 et de 6 telle qu'elle apparaissait precedemment. En 

n 'n 

effet les deux derniers termes de l'expression IX sont identiques 
entre eux selon la condition imposee a la probability de l'amplitude 
du moment magnetique. En consequence les deux equations a deux 
inconnus O n et <J>^ sont solubles. II en resulte que © n et 4> n sont fixes 
par la nature des polynomes en h 1 ] et en g" qu'on a retenus. On peut 
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done caiculer a , etB . en application de l'expression VIII, D'ou 

n-1 r n-l 

on deduit £ , et <j> . . Ainsi de suite on peut caiculer la sequence 

n-1 'n-1 

des diracs qui composent Texcitation. 

A ce stade de Implication on peut resumer la situation en 

5 disant qu'on peut choisir deux polynomes de degre n par mis un 
ensemble important de polynomes possibles (a cause du syteme de 
n+I equations a 2n+2 inconnues). On sait, une fois que de tels 
polynomes ont ete choisis, qu'on est capable de caiculer la sequence 
des excitations sous forme de diracs qui ont conduit, a partir d'un 

10 moment magnetique d 1 orientation donnee, aux polynomes en 
question. En refaisant alors ^approximation par impulsion dure, mais 
a Ten vers cette fois, on peut donner une vaieur pseudo continue ou 
meme continue a Tamplitude de l'excitation. Autrement dit on vient 
de transformer l'expression physique du phenomene. Sur cette 

15 derniere expression on peut effectuer simplement des optimisations. 
Toutes les optimisations porteront sur un choix, en fonction de 
criteres cfoptimisation donnes, des coefficients h 1 ? et g 1 ? dont les 
valeurs sont conformes aux conditions enoncees dans i'expression 
VII, Les calculs cfoptimisation en outre reviennent a *des caicuis de 

20 filtres frequentieis. On sait imposer des conditions cfamplitudes et 
de phases a de tels filtres, 

D f un point de vue rmcroscopique Peffet global de Texcitation 
radiofrequence combinee au mouvement de precession est de 
soumettre Taimantation a une rotation de parametres angulaires © 

25 (w) ; psi (w) et phi (w), ou ©, psi et phi sont les angles cfEuler 
definissant cette rotation, w est la pulsation de resonance 
considered Chacun des trois angles O, psi et phi est fonction, 
comme a et b, de l'historique de l'excitation (t) et peut etre 
trouve, pour chaque w, par integration sur le temps des equations de 

30 Bioch. II existe une correspondance directe entre les angles d ! Euler 
de cette rotation et les amplitudes de probability, a et b, solutions 
de l ! equation de Schrodinger (Albert Messiah, "Mecanique 
quantique", Tome U, pages 918, 919 et 922). Ces correspondances 
s'ecrivent : 



• 1 
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(j phi/2) cos (e/2) exp (j psi/2) 
b = £ exp (-j phi/2) sin (Q) exp (j psi/2) 

avec £ = * 1 



On peut inverser ces relations et obtenir : 



cos & = (aa - bb) = (eta - pp) 



10 soit 



cos 0 = 1-2 = 2 aa - 1 



et 



15 



psi = Arg (a.b) = Arg (z n . a. f>) 
phi = -Arg (a. b) = - Arg (a, p) 



On peut chercher a resoudre le probleme consistant a approximer, 
2Q sur tout un domaine de w, un angle de nutation donne e (w) sans 
s'interesser aux angles psi et phi. 

On sait que si l f excitation est un train de n impulsions de Dirac 
alors les a et£ resultants sont des polynomes d'ordre n en z .11 est 
facile de voir alors, compte tenu du fait que | z | = 1 et done que le 
25 conjugue de z est egal a z l 9 que et aS sont des"polyn6mes" de 
"degre" 2 n (puissance de z de -n a +n). Ces polyndmes ont par 
ailleurs la particulate d'etre symetriques et de pouvoir s'ecrire : 

30 

ou 

Trouver les coefficients Pjqui approximent au mieux une 
reponse en frequence (cos S (w) - l)/2 donnee est un probleme 
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classique de filtrage. Voir par exemple A.V. Oppenheim, R.W. 
Schaefer, Digital Signal Processing, Chapitre 5. On utilisera de 
fa^on preferee l ! algorithme dit cPechange de Remez expose dans 
cette reference et pubiie dans les IEEE Transactions on Circuit 
5 Theory, CT - 19, 189 - 194, Mars 1972 par MM T.W. Parks et J.H. 
McCiellan sous le titre "Chebyshev approximation for non recursive 
digital filters with linear phase". En plus de son caractere cfopti- 
malite cet algorithme presente l'avantage de pouvoir rendre 
compte facilement de la contrainte : 

10 

o^pp<i 

Une fois ies determines on peut trouver les polynomes ^(z). 

En effet si z. est un zero de p alors ^(z) = pp a comme zero z. et 

15 1/z.. A 1' inverse connaissant les 2n zeros de ^(z) on construit j| en 

choisissant parmis chacun des n couples (z., l/zT) un de ses elements 

— li _ 

z^ ou 1/zT. II y a done n choix binaires possibles et, pour un pp 
donne, 2 n solutions possibles pour II en va de merr.e pour a : pour 
aa donne il y a 2 a solutions a possibles. II y a done, pour un cos 0 (w) 

20 approxime au mieux par un ^(z), 4 n solutions pour les 2 autres angles 
cTEuler psi et phi. 

On peut aussi, par recherche des zeros de ^(z) et 1 - £(z), 
examiner toutes les solutions psi et phi, mais on congoit que cette 
operation est longue et fastidieuse. En fait le nombre de solutions 

25 possibles n f est pas aussi grand qu'il y parait de prime abord car 

Talgorithme de Remez tend a placer le plus de zeros possibles sur le 

cercle unite [ z I = 1, et dans ce cas z. - 1/zT. 

I • ■ ii 

En ce qui concerne a on peut generaiement se restreindre a 
chercher une solution a phase minimum. Cette solution correspond a 
30 choisir entre Ies zeros z. et l/zT celui qui est a I'interieur du cercle 
unite. En fait cette solution peut etre trouvee sans avoir a chercher 
les zeros de 1 - f (z) grace a un algorithme simple invente par MM. 
R. Boite et H. Leich ; publie dans Signal Processing, 1981, pages 101 a 
108 par North Holland Publishing Compagny dans un article dont le 
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titre est : "A NEW PROCEDURE FOR THE DESIGN OF HIGH 
ORDER MINIMUM PHASE FIR DIGITAL OR CCD FILTERS". D ! un 
point de vue physique la solution a phase minimum pour a se justifie 
par le fait qu'eile conduit, apres inversion des equations de 
^ ShrodTnger par IX ; a Tenergie de Texcitation la plus faible. En effet 
on sait que parmis toutes les solutions a de aa = 1 -f, la solution a 
phase minimum est celle qui a le coefficient d'ordre zero h^ le plus 
grand. Or en reprenant les relations VI on trouve que : 

n cos n-1 o cos n-1 . cos 

10 h 0 = 2 2 2 

Si Ton prend le logarlthme de cette valeur, et si on considere que 
les contributions apportees par chacun des w. a la nutation totale du 
moment magnetique sont faibles, on peut ecrire, apres avoir 
15 effectue un developpement limite de cosinus © . : 

n 

2 



20 



30 



Log hj - -fE»i) 

" U 



n 



Comme on a vu que etait le plus grand on en deduit que la somme 
des carres des 0. est la plus faible. II en resulte que les impulsions, 
dont la succession provoque une somme des deux angles d'Euler psi + 
phi minimum, sont justement celles qui correspondent a l'energie 
minimale d'excitation pour un angle de nutation B donne. Ce resultat 
n'etait pas connu precedemment . 
25 Une approche plus restrictive de la definition des polynomes a 

et ^ consiste a imposer des le depart des contraintes sur les phases 
arg (a) ou arg (b). En particulier on peut imaginer que arg (b) = 0. 
Ceci est equivalent a imposer que les coefficients du polynome ^ 
soient symetriques, c'est a dire que 



n n 

s . = ^ 

n - i 



Le polynome § peut alors etre trouve par optimisation de Tangle de 
nutation 6 (w) grace a l'alogrithme de Remez. On en deduit alors 
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et aa. On choisit alors a a phase minimum. La figure 2 montre 
Failure d'une impulsion a 180° trouvee en respectant ces conditions. 
Cette figure est associee a un premier tableau, TAB 1, de valeurs de . 
i'amplitude de l'excitation en fonction d'un certain nombre de 

5 periodes (les periodes sont donnees en temps reduit et les amplitudes 
en radian s). Parce que la forme de cette excitation est symetrique 
dans le temps, seules les 1 + n/2 premieres periodes sont tabulees. 
Dans fexemple ii y en a 59. Pour des raisons technologiques cette 
impulsion peut conduire a des puissances cretes, -au centre, trop 

10 fortes pour les generateurs 3 utilises. On peut alors modifier son 
allure en prenant une des autres solutions de a. On aboutit alors a 
une impulsion dont Tenveloppe est celie representee sur la figure 3, 
et qui correspond t dans les memes conditions, a un deuxieme tableau 
de valeurs : TAB 2. A Texperience il apparait que le respect a 3 ou 5 

15 % pres des valeurs indiquees dans ces tableaux rfest pas prejudi- 
ciable au resultat des experimentations de RMN qui les utilisent. 
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TAB 1 



TAB 2 



-3. 90491O8882Oe2lE-OO2 
- 6 . *663 1 9801 34976E -003 
-9. 231 66598S98407E-003 
1 96774 630407086-003 
- 6 . 8 34 328 25998 1 23E-003 
-6. 12591456552324E-003 
• 7 . 0596 1 652094723E-003 
-3. 63477537258307E-003 
-3. 8625391 1015924E-003 
-1 .7e4664 35l3e458E-003 
6. 25476838024989E-OO4 
"! 2644069Ol4«6B^E-OO3 
6 . 095830 76591 332E-0C3 
9 . 0521 6904390593E-O03 
1 . 20555 1 4256 1 780E-002 
l.Wi 967 68 191 B74E-002 
1 . 783 1 74 1 3596757E-002 
2 . 04537546376239E-002 
2. 272461 68283321E-002 
2 . 456301 46956734E-002 
2 . 58 7 1 028 7 1 5364 2E -002 
1 . 6557 866 1 205690E -002 
2.6543135565A137E-002 
2. 57591 1 39192740E-002 
2.41 545403 ! 35 1 47E-002 
2. 1 69569 17805766E-O02 
1 • 837 1 3380832749E-OO2 

1 . 41 945062 109303E-002 
9. 2045766O885527E-rOO3 
3. 4676S560648O96E-0O3 

. - 2 . 92329 1 00733593E-003 
-9. 8473390S633224E-O03 
-1. 715301 73784658E-002 
-2. 4637153581 6664E-002 
•3. 2 144221 94642 14E-002 
-3. 9369381 401 6602E-002 
'-4 . e062 1 260029 1 37E-OC2 
-5. 19305 174B98276E-002 
-5. 663646967486 10E-OO2 
-5. 98922285225749E-002 
-6. 13121319539390E-002 
•6.0598114077B251E-002 ' 
-5- 744947391 12293E-002 
•5. 15620362699079E-O02 
-4. 26629999764663E-002 
-3. 0305568 1 138646E-002 
- 1 . 486 1 8 1 0837340 1 E-002 
4. 48331 83831 7*62E-003 

2. 773 1 4 167977363E-002 

3. 30461 46234 171 3E-002 
8. 64972903264657E-002 
K21931S61422393E-001 
1 . 60897063098379E-O01 
2. 0220306397 1 342E-O01 

2. 4385907742 1071E-001 
2.82823939221 31E r 001 

3. 151 869792 18870E-001 
3. 3682B7B0B39443E-OO1 
3. 445 19 lOl 220846E-001 



-6. 43620 15321 999 IE -002 

• 2. 229 067 39796626E-002 
-2. 49869020943 108E-002 
•2. 68990157438880E-002 
-2. 77707368058623E-0O2 
-2. 7 3394 7 1 1 9 1 63 1 1 E -002 
-2. 54479434557955E-002 
-2. 18635472488279E-002 
-1 . 64935463610009E-002 
-9.2973U66463141E-003 
-3. 19132207764246E-004 

1 .030290774 08931 E-002 
2. 2324 9500506 28 1 E-002 

3. 539O662723737BE-O02 

4 . 90274 1 906 1551 3E-002 
6.26451 14B417868E-002 
7 . 5346637 1 96 7290E -002 
8. 69479228131 243E-002 
9.601 5522 1 4 1 5 1 1 6E-002 
1 .01921802249669E-001 
1 . 039 1 3 1 384 39 1 67E-00 1 
1.01 380 1878897788-001 
9. 39 13552 1 524 1 7 BE -002 
B. 13463363002227E-002 
6 * 3739022 ! 677407E-^02 
4. 13790603273033E-002 • 
1 . 4 74863481 49 16 1 E-002 

• 1 . 54 839234400 1 07E -002 
-4. B44244 1 01 6884 1E-002 
-B. 2985233 168 4 B52E -002 
-1.1 7623793321 197E-O01 
- 1 . 3047 1 463648678E-001 
-1 . 79244 B67442907E -001 
-2. 01 439013842493E-001 
-2. 1 47692836 13060E-001 
-2. 17434784481 U7E-O01 
-2 . 087 1 69 1 8900338E-00 1 
-1 . 8903007 7 06634 5E -001 
- 1 . 597389703706 1 6E-001 
- 1 . 226474076B6439E-001 
-8. 04 03761 297993 1E-002 
-3. 4389B606033687E-002 

1 • 34844833846245E-002 
6. 167 4285251 404 1E-002 
1 . 085826O6007642E-001 

1 • 3237236040^07 1E-001 

1 . 9OB968709850O5E-OO1 
2.218411 1 22204 1 9E-00 1 

2. 431 14767635960E-O01 
2. 33390364 178345E-O01 
2.52376B07632722E-001 
2.41 983370764 297E-O01 

2 • 24 0637053 1 1 999E-00 1 
2.01994161 175123E-001 
1.7 905167S99749E-O01 
1 . 37651 976499981 E-001 
1 .40485395387081E-001 
1 . 29003722671 106E-001 
1 . 24 1 64 37623286 1 E-001 
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REVENDICATIONS 

1 - Procede cfexcitation radiofrequence dans une experi- 
mentation de RMN, excitation en reponse a laqueile les bascu- 
lements de moments magnetiques de particules soumises a i'experi- 
mentation sont definis en amplitude et en phase pour un ensemble de 
frequences de resonance donne, caracterise en ce que s 

- on ecrit ^equation de Schrodinger relative a I'amplitude de 
probabilite de la fonction cPonde de spins resonnant a une frequence 
de resonance donnee et dont i'hamiltonien depend de la frequence et 
de l'excitation radiofrequence qu'ils subissent, 

- on integre cette equation par rapport au temps, 

- on en deduit les expressions analytiques polynomiales des ampli- 
tudes de probabilite a la fin de l'excitation en fonction de la 
frequence de resonance, ces expressions polynomiales etant de degre 
egai a un no rubra de periodes elementaires dont la somme des durees 
est egale a la duree de l'excitation, et ces expressions etant 
parametrees par des coefficients representatifs de l'historique de 
1 ! excitation pendant ces periodes, 

- on optimise ces coefficients a l'aide de Talgorithmes de calcui de 
fagon a ce que les polynomes s'approchent au mieux rfune forme 
desiree, 

- on deduit les caracteristiques pendant chaque periode de Texci- 
tation qui conduit a cette reponse definie. 

- et on excite les particules avec une telle excitation, 

2 - Procede selon la revendication 1, caracterise en ce que 
22 pour integrer on eff ectue une approximation par impulsions dures, 

3 - Procede selon l'une quelconque des reyendications 1 ou 2, 
caracterise en ce que pour optimiser on recherche les expressions 
polynomiales a dephasage minimal qui correspondent a une exci- 
tation la moins puissante possible. 

4 - Procede selon la revendication 3, caracterise en ce que 
pour exciter on exite avec une une excitation dont Failure est celle 
de la figure 2 representative du tableau 1. 
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5 - Procede selon la revendication 3 caracterise en ce que 
pour exciter on excite avec une excitation dont l'ailure est celle de 
la figure 3 representative du tableau 2. 
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